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Nowadays, geometric tools are being used to treat a huge class of problems of quantum information science. 
By understanding the interplay between the geometry of the state space and information-theoretic quantities, it 
is possible to obtain less trivial and more robust physical constraints on quantum systems. Here we establish a 
geometric lower bound for the Wigner-Yanase skew information (WYSI), a well-known information theoretic 
quantity recently recognized as a proper quantum coherence measure. In the case of a mixed state evolving 
under unitary dynamics generated by a given observable, the WYSI between the state and the observable is 
bounded from below by the rate of change of the state’s statistical distinguishability from its initial value. Our 
result shows that, since WYSI fits in the class of Petz’s metrics, this lower bound is the change rate of its 
respective geodesic distance on quantum state space. The geometric approach is advantageous because it raises 
several physical interpretations of this inequality under the same theoretical umbrella. 
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I. INTRODUCTION 

Coherence is a striking feature of the quantum realm due to 
interference phenomena [1], In fact, it is in equal footing with 
entanglement and other correlations whose meaning evades 
the classical view. Although quantum optics has proved to 
be a fruitful branch for quantum coherence studies [2], recent 
results suggest its connection with thermodynamics [3] and 
quantum biology [4], Even some condensed matter phases, 
such as superconductivity and its emergent properties, dis¬ 
play signatures of quantum coherence [5]. As we know, the 
inability to perform basic tasks in quantum information pro¬ 
cessing is often related to coherence loss and, because of that, 
the interplay between noise and decoherence still holds as a 
key challenge in open quantum systems. Recently, Bromley 
and coworkers [6] reported a possible way to circumvent this 
problem. Summarizing, they found a regime called freezing 
conditions in which coherence remains unchanged during the 
nonunitary dynamics. 

Despite its fundamental role in physics, there is no unified 
way to characterize and quantify coherence. In consonance 
with results presented in Ref. [7], a recent approach due to 
Baumgratz et al [8] established a new paradigm in this sce¬ 
nario. By employing a rigorous mathematical framework for 
identifying proper coherence measures, they were able to clas¬ 
sify natural candidates for coherence quantifiers based on dis¬ 
tance measures, particularly relative entropy, l p -norms, and 
fidelity. Simultaneously, Girolami [9] proposed another quan¬ 
tum coherence measure based on the Wigner-Yanase skew in¬ 
formation (WYSI), which shares the same reliable criteria of 
Ref. [8]. Besides the theoretical background, this work of¬ 
fers an efficient route to experimentally access the quantum 
coherence of an unknown state. 

In the present work, we focus on skew information to pro¬ 
vide an information-geometric lower bound for coherence 
measures. Introduced by Wigner and Yanase half century 
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ago [10], skew information 

J(p,70:=-iTr([^,7C] 2 ) 0) 

is a measure of the non-commutativity between a state p and 
an observable r K. Operationally, this quantity is deeply more 
interesting than other coherence quantifiers because its cal¬ 
culation does not involve any optimization techniques. Also, 
it describes a constant of motion in closed quantum dynam¬ 
ics when the observable 'K is a conserved quantity, i.e., it 
commutes with the hamiltonian generating the evolution of 
the system [11], Furthermore, WYSI is nonnegative, con¬ 
vex and vanishes if and only if the state and the observable 
commute [12], It is also bounded by the variance of C K, 
I(p,'K) < ('K 1 ) p - ( r K)p, an interesting property discovered 
by Luo which also noticed that the inequality is saturated for 
pure states [13]. This measure was later generalized by Dyson 
as 

I p (p,'K)-.= - l -Tx{[p p ,<K][p l - p ,‘K)) (2) 

with 0 < p < 1, being called Wigner-Yanase-Dyson skew 
information (WYDSI), and its convexity proved by Lieb [14], 

There are several interpretations of the skew information, 
each one related to a particular viewpoint of the quantum be¬ 
havior. Actually, the original one discusses the uncertainty in 
the measure of observables not commuting with a conserved 
quantity - basically, the content of Wigner-Yanase-Araki the¬ 
orem [15]. Similarly, WYSI supports a new type of Heisen¬ 
berg uncertainty relation [16], quantifies the quantum uncer¬ 
tainty of local observables [ 1 7] and has applications in quan¬ 
tum reference frames and metrology [18]. It is also possible to 
detect entangled states through a Bell-type inequality derived 
from the skew information [19]. 

WYSI is also an asymmetry measure, i.e., it quantifies sym¬ 
metry breaking in a given state [20], This is a promising sub¬ 
ject in quantum information which finds support on the asym¬ 
metry theory and classifies coherence as a resource [21]. In 
this context, Noether’s theorem is a powerful tool to character¬ 
ize conservation laws from symmetries in closed quantum sys- 
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terns because each asymmetry measure is a conserved quan¬ 
tity. Nevertheless, recent efforts have elucidated some asym¬ 
metry properties of pure states and quantum channels [7], but 
the mixed state case is rather complex and less exploited. This 
happens because, when dealing with mixed states one must 
search for conservation laws which are not captured in its 
essence by Noether’s theorem [22]. As advocated by Marvian 
and Spekkens [23], an asymmetry measure based on WYSI 
could fill this gap providing a way to point out more subtle 
features of conserved quantities. 

The main result of our work is that, for closed quantum sys¬ 
tems, the skew information, J(/y,A^), between an evolved 
mixed state, p^, and the observable, K^, generating its evo¬ 
lution is lower bounded by the rate of change of the distin- 
guishability between the evolved and the initial, po, mixed 
states 


— cosUXpo.Py,)] 

dip 
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Here -CipQjpf) is the Hellinger angle between the initial state 
and the evolved state. The evolution is given by a family of 
unitary transformations U v which changes continuously with 
respect to the parameter < p. The observable may or not 
depend on the parameter <p and it is connected with the op¬ 
erator Uy through the relation /Q = -ihlj^idul/dip). In our 
approach, the encoded parameter tp can assume different in¬ 
terpretations depending on a specific physical situation. For 
instance, it could be the phase difference introduced in an in¬ 
terferometric protocol, with K v being the generator of the ro¬ 
tation, or the time in a dynamical evolution, in which case K v 
would be the Hamiltonian of the system. 

The paper is organized as follows. In Sec. II we review 
the necessary and sufficient conditions that WYSI should sat¬ 
isfy in order to be a proper coherence measure. In Sec. Ill 
we point out that this information-theoretic quantity defines a 
monotone Riemannian metric due to the Petz’s theorem [24] 
and its respective geodesic distance on quantum space state 
is given by the Hellinger angle [25]. By exploring these two 
interfaces, in Sec. IV we demonstrate an inequality which as¬ 
signs a geometric meaning to coherence measures. In other 
words, we show that, for closed quantum systems, the skew 
information between an evolved state and the observable gen¬ 
erating the evolution is lower-bounded by the rate of change 
of the distinguishability between the evolved and the initial 
states of the system. In Sec. V we provide an example in or¬ 
der to illustrate our claim. Finally, in Sec. VI, we present our 
conclusions. 


n. WYSI AND QUANTUM COHERENCE MEASURES 


notions, in the following we will present these ideas in a rig¬ 
orous fashion based on Refs. [8, 9]. 

It is well known that quantum operations are described by 
dynamical maps, i.e., quantum channels. In particular, the ac¬ 
tion of completely positive and trace preserving (CPTP) map 
& on the state p can be synthesized as &(p) = Yj^K M pK^, 
where j K jt } is a set of Kraus operators satisfying A ; ', K jt = I. 
Let 7/ be a finite-dimensional Hilbert space with d — dim 1 //. 

Choosing a fixed basis [|i)]i=i.</, the subset of incoherent 

states J? c 77 encompasses those whose density matrix is di¬ 
agonal in this basis. So, an incoherent channel (ICPTP) is the 
map c ,A for all p, i.e., transform incoherent states 

into incoherent states. In other words, this constraint excludes 
any coherence generation process. 

As demonstrated by Girolami [9], skew information is a 
faithful coherence measure since it satisfies the axiomatic pos¬ 
tulates proposed by Baumgratz et al [8], First, it is convex, 
non-negative and vanishes for all incoherent states p e J?. In¬ 
deed, I(p, < K) — 0 if and only if [p,7C] = 0, i.e., state and 
observable can be diagonalized simultaneously. Secondly, it 
is monotonically nonincreasing under ICPTP maps and does 
not increase on average under a von Neumann measurement, 
J(p,7C) > ^PiuIiKiipK^/K), where /?,, = Tr (K^pK^). 


III. WYSI, PETZ METRIC AND HELLINGER ANGLE 

WYSI is a robust information-theoretic quantifier due its 
enormous versatility. Actually, skew information also can 
be interpreted from a geometric perspective. The most re¬ 
markable approach to achieve this goal is indubitably due to 
Morozova-Cencov [26] and Petz [24, 27], by using mono¬ 
tone metrics on the quantum state space. In this space the 
set of density operators (p > 0 and Trp =1) constitute a dif¬ 
ferentiable manifold equipped with a suitable monotone Rie¬ 
mannian metric. By monotone metrics we consider the ones 
that are defined by positive, continuous and sesquilinear inner 
products which are also contractive under CPTP maps. 

The Morozova-Cencov-Petz theorem provides a friendly 
way to demonstrate that skew information fits in the category 
of monotone metrics and describes a particular kind of quan¬ 
tum Fisher information [28]. Generally speaking, the theo¬ 
rem states that there exists a bijective correspondence between 
monotone metrics and operator monotone functions given by 

gf(A, B) := Tr[A c/(Jf ,3£)E], (4) 

where A and B are traceless hermitian operators and 

c f(.x>y) ; = J , . ( 5 ) 

yf(x/y ) 


In order to characterize skew information as a coherence 
measure, it is essential to establish the concept of incoherent 
states and incoherent operations. An incoherent state is one 
that has no coherence, i.e., its off-diagonal elements are equal 
to zero. In the same way, an incoherent operation is one that 
does not create any kind of coherence. Despite the intuitive 


is a symmetric function, c/(x,y ) = c/(y,x), and fulfils 
Cf(ax,ay) = a~ l Cf(x,y) with x,y > 0. Here PAG = pG 
and PAG = Gp are commuting operators, [Jzf, PA\ G = 0. 
Besides, the function fit) is if) operator monotone, i.e., for 
any density matrices ,e/, PA such that 0 < .&/ < PA, then 
0 < f(sf) < fi-P$)\ iii) self-inversive, /(f) = f/(l/f); and iiii) 
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normalized, /(1) = 1. Naturally, there are as many mono¬ 
tone metrics as there are operator monotone functions, which, 
according to Petz, represent a vast garden of monotone met¬ 
rics [27]. 

As pointed out by Gibilisco and Isola [29], an ordinary ele¬ 
ment of the tangent space of the density matrices manifold is 
given by i[p,'K\, where 7C is an Hermitean operator. Notably, 
choosing A - B - i\p, < K\ and taking /(f) = (1/4)(V?+ l) 2 
such that Cf(x,y ) = 4/( y/x + r/y ) 2 , it follows that 

g f (i[p, <K), i[p, 7C]) = 4Tr (i[p, <K] Jt? i[p, 9Q) 

= 4{JZ / 1 [p, 7f ], ,M~ 1 [p, 7f] >, (6) 

where ± y/W and (A, B) = Tr(A 1 B) is the 

Hilbert-Schmidt inner product. Since [p, 7C] = (22? - = 

. # + ./# ( K, we get the monotone metric 

g f (i[p, 7C], i[p, 7f]) = JfJK) 

= 4 <[yfi,9C],[y/p,X}> 

= -4Tr ([^,7C] 2 ) 

= 8 J(p, 70, (7) 

which, up to a constant factor, is exactly the Wigner-Yanase 
skew information. Recalling the multiple facets that WYSI 
embodies, Eq. (7) indicates a clearly connection between 
coherence measures and information geometry. It is worth 
mentioning that other authors also addressed the geometri¬ 
cal features of WYSI in a rigorous viewpoint [30]. Recently, 
Brody [31] has demonstrated that the space of pure and mixed 
states is equipped with a dual metric structure which assign a 
clear meaning to the WYSI in the geometric realm. Besides, 
his approach also enabled to derive corrections to the Heisen¬ 
berg uncertainty relation based on skew information. 

Since the quantum state space is endowed with a metric 
structure, it is natural to ask about distances, curvature and 
other geometric properties. Particularly, the notion of distance 
between states has been the subject of discussions initiated 
decades ago under the spotlight of statistical inference [32], 
In a pioneering work, Wootters employed the statistical dis¬ 
tance concept as a proper distinguishability measure between 
statistical probabilities [33]. The geometrization of this prob¬ 
lem emerged years later with Braunstein and Caves [34] who 
defined a Riemannian metric and its respective line element 
ds from a suitable distinguishability quantifier between close 
states. Though their description was based on a physical 
ground, it is analogous to that one developed by Petz which 
relies on monotone metrics. Summarizing, the main message 
about those works lies on the close relation between state dis¬ 
crimination and geometric distances. 

Following Petz’s approach for the WYSI monotone met¬ 
ric, it has been proved that the distance between two den¬ 
sity matrices on quantum state space is D{p,cr) — 2 - 
2Tr( yfp a/o") [25]. This quantity is the quantum analogue 
of the classical Hellinger distance [32]. Our discussion on 
the geometric properties of WYSI should include a few lines 
about geodesics - the shortest distance between two den¬ 
sity matrices on the quantum state space - associated to the 


Wigner-Yanase monotone metric. It was shown that the cor¬ 
responding geodesic distance joining the density operators p 
and cr is given by the Hellinger angle [25, 35] 

-C(p, cr) = arccos[Tr( yfp y/cr)] . (8) 

The quantity JH(p, cr) = Tr( yfp y/cr ) is called quantum affinity 
and describes how close two states are on the quantum state 
space [32], Moreover, it is remarkable that quantum affin¬ 
ity is bounded from below by the Quantum Chernoff Bound 
(QCB) [36], 


IV. GEOMETRIC LOWER BOUND ON QUANTUM 
COHERENCE 


We now provide a lower bound for the quantum coher¬ 
ence measure based on the skew information. Let us fo¬ 
cus on a driven closed quantum system described initially by 
a mixed state po which undergoes a unitary transformation 
Pip = U vPo ul. Essentially, this operation encodes the param¬ 
eter ip on the input state and does not change its purity. The 
operator LQ characterizes a family of unitary transformations 
labelled by <p. Besides, it is worth to mention that U v changes 
continuously with respect to this parameter. The reason for 
starting from a mixed state is twofold: first, the skew infor¬ 
mation is bounded by the variance when dealing with mixed 
states [11, 12], Actually, this result was improved later by a 
variance lower bound which is tighter than this one based on 
the skew information [31], Moreover, it also allowed to de¬ 
rive an entire family of higher-order corrections to the uncer¬ 
tainty relation supported by WYSI by exploiting its connec¬ 
tion with the quantum analogue of the conditional variance; 
second, because all Petz’s metrics - particularly the Wigner- 
Yanase one - becomes the well known Fubini-Study metric 
for pure states [37]. 

Considering the Wigner-Yanase metric in the quantum state 
space, according to Eq. (8) we obtain 


— cosUXptnP^.)] 
dip 


—Tr( Vpo yfp/) 


(9) 


Since U v is a unitary operator, it is possible to write yfp/ = 
Utp y/poU/p (see Appendix A) which implicates the von Neu¬ 
mann equation 


d -dp/ i 

dO) 

where we used that ( dU lf ldip)U] f — — U/{du/p/dip) and defined 
the Hermitian operator 


K ip = -ihU l p—^ . (11) 

dip 

In general, the operator depends on the parameter < p. How¬ 
ever, it is worth to notice that as a special case, when the ob¬ 
servable Kp is independent of ip, i.e., - K, thus the unitary 

evolution is given by — e~ llfK . 







4 


Returning to the general case and substituting Eq. (10) into 
Eq. (9), we have 


d 

dip 


COsUXpo.Pp)] 


j: |Tr( Vpo[^> A^])| • 


( 12 ) 


Equation (12) is the starting point for establishing the lower 
bound on the skew information. Actually, this goal is reached 
by noting that 

|Tr( yfp&K,, ^])| < II yfmiWlK,, ylp^h . (13) 


where we used the Cauchy-Schwarz inequality |Tr(A7?)[ < 
||A|| 2 || 5 || 2 , with ||A ||2 = y/TriA^A) being the Schatten 2 -norm 
(also known as Hilbert Schmidt or Frobenius norm) [38], 
Combining Eq. (13) with || ^polb = 1 and substituting the 
result into Eq. (12), we obtain 


— cos[£(po,Pp)] 

dip 




(14) 


On the other hand, note that 


\\{K V , yfp^Wi = ^-Tr([ y/p^, K\p] 2 ) = ^2lip^K ^), (15) 


where I(p^,K£> = —(l/2)Tr([ y[p^, K v ] 2 ) is the Wigner- 
Yanase skew information between the evolved state p v and the 
observable K^ that generates the dynamics. Therefore, substi¬ 
tuting Eq. (15) into Eq. ( 14) we obtain a lower bound in terms 
of WYSI and Hellinger angle as follows 


— cosUXpo.p^,)] 
dip 



(16) 


Eq. (16) is the main result of this paper. It is important 
to highlight that it encompasses any class of continuous uni¬ 
tary transformations U^ indexed by the parameter ip, as well 
as initial and evolved mixed states. As a particular case, 
recalling that a unitary evolution does not change the pu¬ 
rity of a quantum state, if po is pure, then will also be, 
and the skew information reduces to the variance of K^, i.e., 
Iiptp, K v ) = (A Ky,) 2 = (X -2 ) - {K tf ) 1 . In this regime, the lower 
bound becomes 


A K v 




(17) 


where f(ip) = Tr[pop^]/Trp 2 defines the relative purity, which 
played a special role for the investigation of quantum speed 
limits under the closed dynamics [39]. 


V. EXAMPLE 

To illustrate the use of the bound indicated in Eq. (16), we 
now consider the single qubit case. Let po = (l/2)(7 + r t) ■ <f ) 
be the initial state (7 denotes the 2x2 identity matrix, ?o is a 
3-dimensional vector with |r () | 2 = < 1 and & = {o - !, <x 2 , cr 3 ] 

is the vector of the Pauli matrices). The dynamics is governed 
by the self-commuting local observable K v = mini + • r?j, 


i.e., [Kp, K^] = 0 for all < p and ip', where m and a are posi¬ 
tive constants and is an unit vector, In^l = 1. The system 
evolves under a general unitary operator U v given by 


Uy = exp 



= e ,Sa [I cos y - i(t,,p ■ (?) sin 7 ] , 


(18) 


where 7 = mtp\I. i p\/fr and 6 - zuip/h are dimensionless con¬ 
stants and also 



Essentially, the initial state po undergoes the unitary trans¬ 
formation p, f — UpPoUl which encodes the parameter ip. 
It is possible to verify that the final state can be written as 
p<p = (1 /2)(7 + • (?), where (see Appendix B) 

= cos(2y)r 0 + [1 - cos(2y)](I^ ■ r 0 )i^ + sin(2y)(I^ X r 0 ) . 

( 20 ) 

The vector keeps whole information about the parameter p 
and has the same magnitude as the initial vector ?o, i.e., \r^\ = 
|r 0 l - '"(). Particularly, as a special case, if ti v is independent of 
the parameter^, i.e., = h, then Eq. (19) implies that = h 

and, consequently, 7 = 6 — mp/ti. 

In order to calculate the Hellinger angle we need to de¬ 
termine the trace of the product of operators yfpo and \jjp- 
The analytical expressions for the square root of a single qubit 
state can be found in Appendix C. Since the modulus of Bloch 
sphere radius remains constant under the unitary transforma¬ 
tion, it is possible to verify that the cosine of the Hellinger 
angle becomes 

cos[-£(p 0 ,p^)] = + £-(r<p • r 0 )] , ( 21 ) 

where <f ± = 1 ± 1 - r 2 is independent of the param¬ 

eter ip. From this result is straightforward to check that 
dcos[£(p 0 ,p<p)]/dip - (g-./2)[(dr v ,/dtp) ■ r 0 ]. Similarly, the 
Wigner-Yanase skew information is given by 

I{p v , K^) = nr 2 ^_|^ X n ^| 2 . (22) 

Substituting the derivative of Eq. (21) on the parameter ip 
and Eq. (22) into Eq. (16), we finally obtain the bound 
■\fj[l\{dr< P ldip) - Pol <2 V2(nr/?z)|f^ x h^\. To clarify, choos¬ 
ing the probe state po = (1 - ro)7 + ro|^)(^|, | if/) = 
(1/V2)(|0) + 11 )), and h = (0,0, 1 ), which corresponds to 
take ?o = ro(cos (f>, sin <p, 0 ) (0 < ro < 1 and 0 < <p < 2n) 
and = m(al + cr z ), the cosine of the Hellinger angle is 

cos[.£(po,Py,)] = (l/2)[£+ + £_cos(2ip)] and the WYSI gives 
J(p^,7f^) = nr 2 £_sin 2 0. Combining both results we obtain 
the bound V^I[ sin(2ip)| < V2(ur/^)| sin 0|. It is interesting 
to note that, while JXpo>P<p) is a function of the parameter ip, 
WYSI is independent of this phase and describes a constant 
of motion during the unitary evolution. It is worth mentioning 
that although we focused attention in the single qubit case, our 
calculations can be extended to a system of N qubits. 
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VI. CONCLUSIONS 
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In this article we established a geometric lower bound for a 
proper quantum coherence measure based on Wigner-Yanase 
skew information. This information-theoretic quantity is re¬ 
garded as a particular extension of Fisher information and 
can be seen as a monotone Riemannian metric due to the 
Petz’s theorem [24], Moreover, its related geodesic distance 
on quantum space state is given by the Hellinger angle [25]. 
In opposition to many other distance measures such as Bu¬ 
res angle or even relative entropy, Hellinger angle is advan¬ 
tageous quantity because is technically easier to calculate and 
more intuitive to obtain from its classical statistical analogous. 
Despite those motivational issues, it has received little atten¬ 
tion beyond that devoted to the exploration of its useful al¬ 
gebraic properties to the information theory. It is important 
to emphasize that our result shows that, since geodesic dis¬ 
tance quantifies the discrimination of two density operators in 
the context of quantum statistical estimation theory [32], skew 
information is bounded from below by the rate of change of 
distinguishability between two states on quantum state space. 

Our result opens a wide range of possible physical inter¬ 
pretations. First, inequality Eq. (16) suggests a route for bet¬ 
ter understanding the phase estimation paradigm in quantum 
metrology [40]. In fact, it can provide a precision bound for an 
unknown parameter <p encoded by the unitary transformation 
in the initial state. Therefore, the bound essentially depends 
on the derivative of the Hellinger angle with respect to this 
parameter. 

In particular, choosing ip - t, where r is time, it can 
be shown that our inequality gives rise to a new quantum 
speed limit [41], In contrast with the original one proposed 
by Mandelstamm-Tamm [42], and later generalizations for 
driven closed systems [43], this speed limit depends on WYSI 
and the Hellinger angle rather than the Bures angle or the vari¬ 
ance of the hamiltonian. 

Besides, it seems possible to attribute a thermodynamic 
meaning for this bound by investigating the connection be¬ 
tween nonequilibrium entropy production [44] and the ther¬ 
modynamic length [45] involving quantum protocols at finite 
temperature. This could provide a thermodynamic interpreta¬ 
tion for the existence of the quantum speed limit. 

Finally, in a future work it will be crucial to investigate the 
eventual relation between geometric bounds and the univer¬ 
sality class of Petz metrics which fulfils the requirements for a 
quantum coherence measure. Moreover, to enlarge the present 
analysis, take into account the open quantum dynamics would 
be essential not only for the foundations of quantum infor¬ 
mation theory but also for realizing quantum technology in a 
noisy scenario. From the experimental point of view, by ex¬ 
tending our conclusions to N quibt systems, the bound in Eq. 
(16) could be experimentally investigated through a measure¬ 
ment scheme based on two-point correlation functions [46]. 
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APPENDIX 

A. MATRIX POWERS 


In this section we will demonstrate that the identity 
(VAVQ' S = FA'F 1 still holds for 0 < s < 1, where A is a pos¬ 
itive matrix (A > 0) and V is an unitary operator, V 1 = V 1 . 
In order to reach the main goal, let us consider a monotone 
function f(a) = a s for a > 0. It can be demonstrated that f(a ) 
has the following integral representation [38] 


a 


S 


sin(7Ti) 
7 TS 



a 

a + x 


dp(x ), 


(A.l) 


where p(x) - X s , dp(x) = sx s ~ l dx, is a positive measure on 
(0, oo). This relation can be extended to the positive and non¬ 
singular operator A as follows [47] 

A s = Sm(nS) f A(A + xlf'dnix) . (A.l) 
ns Jo 


Actually, the last condition can be relaxed if A is full rank or 
if we assume that the inverse operation is taken on its suport, 
i.e., the vector subspace spanned by the eigenstates with non¬ 
zero eigenvalues [48]. Considering the transformation V / AV /t , 
follows 


(VAF f + xiy 1 = [V(A + Jr /)V t ]“ 1 

= V(A + x iyW f . (A.3) 


Therefore, we can prove our main goal combining the pre¬ 
vious equality with the integral representation indicated in 
Eq. (A.2), i.e., 

smins) r*°° 

(kAk 1 )’ = —— VAV\VAV f + xiy'dp(x) 

ns Jo 

= V Sln(;r ‘ y) f A(A + xI)~ l dp(x)V f 
ns Jo 

= VA s V f . (A.4) 


Particularly, given the evolved state p^ = U^poU^, choosing 
A = po as the initial mixed state and the unitary operator V = 
U v , for ,s=l/2 Eq. (A.4) allows to demonstrate the relation 

x/p<p (j ip y~poU<p . 


B. UNITARY EVOLUTION 

In this section we describe the calculation of the evolved 
state p v in the single qubit context. Let us assume that the 
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quantum system dynamics is governed by the local observable 
K v = m(al + hy ■ <?), where m and a are positive constants 
and hif is an unit vector, i.e., \h v \ - 1. By hypothesis, this 
observable is self-commuting, i.e., [ K\ p , K^] = 0 for all tp and 
ip'. The system evolves under a general unitary operator 
given by 


U, f = exp 


i r 


e l6a exp [-iy(t v ■ (?)] 

'*''[/ cos y - ift.y ■ (?) sin y] , 


: e 


(B.5) 


where y = vj^lip/h and 6 = zmp/Ti are dimensionless con¬ 
stants and also 



In particular, if h^ is independent of the parameter < p, i.e., h lf — 
h, then %p = h. Returning to the general case, let be an initial 
single qubit mixed state po = (1 / 2 )(/+ ?o • c?), where I denotes 
the 2 x 2 identity matrix, |?o | 2 = r^ < 1 and & = {cri, <X 2 , cr 3 ) 
is a vector of the Pauli matrices. The probe state po undergoes 
the unitary transformation p^ = U if po U]p and can be written as 


P v = 2 ^ + <To • <?)cos 2 y + i[r 0 ■ (?, l v ■ (?] sin y cos y+ 

+ ■ (?)(? 0 • <?)(£* ■ cr)sin 2 y) . (B.7) 


and (Aif)ji '■= €jki(£<p)k- The matrix element (A v )j[ satisfy the 
identity 

(A v )p(A ¥ ,) s ; — € s jk€ s ^i(Xip)k(Xip)^ 

= (SjfjSki - 8ji6kii)(% ip )k(J‘ip)n 

= (% P ) J <X v )i-S J i, (B. 12 ) 

where we used the Einstein summation convention and the 
property (XflkXftk = |2 ^| 2 = 1. From this expression we 
have (£ v ) j(X<p)i = 6ji + Aj S (tp)A s i(ip) and therefore 


(S ip)ji = Sji + [1 - cos(2y)](A 40 ) JJ (A v ) J ; + sin(2y)(A^) j7 . 

(B.13) 

The matrix A^ is called skew tri-idempotent because fulfils 
A 2 = —A ip . This property can be verified starting from the 
triple product 


(A,)p(A,) 

s/d^^-ip)/ul — id(^(p) a VJj 

= “(Ay) j, . (B.14) 


Note that the last equality in the expression above was ob¬ 
tained by using the identity Z • ( 2 ^ x 2 ^) = £i, m X^)nX v )n = 0 . 
From the result obtained in Eq. (B.14) the matrix S v can be 
written as 


S,p - I + [1 - cos(2y)]A 2 + sin(2y)A^ 
= e 2 rA . 


(B.15) 


Exploring the algebraic properties of Pauli matrices is possi¬ 
ble to check that (a ■ f)(b ■ (?) = (a- b)I + i{axb)■ tf. Combin¬ 
ing this relation with the vector identities a ■ (a xb) = 0 and 
ax(b xc) = (a ■ c)b - (a ■ b)c, we obtain 

[?o • (?, 2 v • (?] = -2iXip X r 0 ) • (? (B. 8 ) 

and 


From this relation is possible to identify the explicity form of 
matrix A v . First, this matrix has all diagonal elements equal 
to zero, i.e., (A v )jj = ejkjXftk = 0. Second, the matrix A^ 
is anti-symmetric because (A^) /; = ei k] (X<f )k = SjkiXftk = 
-(A v )ji. On the other hand, given that (A ^) 23 = 6213(^)1 = 
-(£?)i> (A ¥ ,)i 3 = 6123(^)2 = ~Xf)i andlA^)^ = ei^Xiph = 
-( 2 ^) 3 , is immediate to write A^ as 


tip ■ (?)(?o ■ &)(%, ■ <?) = [2(tp ■ ?o)Xp - ? 0 ] • c? ■ (B.9) 


A <p = -iX ip -f, (B.16) 


Substituting Eq. (B.8)-(B.9) into Eq. (B.7) and performing 
the calculations, it is possible to verify that the evolved state 
becomes p v = (l/2)(7 + r,■ (?), with 

r,' f = cos(2 y)r 0 + [1 - cos(2y)](l; v ■ r 0 )2^ + sin(2y)(l:^ X r 0 ) . 

(B.10) 

It is worth to emphasize that both vectors p, and ?o has the 
same absolute value, |?^| 2 = [?o | 2 = r In other words, the uni¬ 
tary transformation does not change the modulus of Bloch 
sphere radius during the dynamics. 

Now we will provide another proof for Eq. (B.10) which 
is based on the Rodrigues’ rotation formula. Summarizing, 
through this approach the vector /y is completely determined 
by the action of a rotation matrix on the initial vector r ( ). In 
order to understand this property, let be the j—th component 
(r<p)j = (S ji(r 0 ),, with 


where J = is a vector whose components 

are given by the generators of the adjoint representation 
(3-dimensional) of SU(2) algebra, 


0 0 0 


0 0 1 


0 -i 0 

0 0 -i 

> Jl = 

0 0 0 

. ?3 - 

i 0 0 

0 i 0 


-i 0 0 


0 0 0 


(B.17) 

Finally, substituting Eq. (B.16) into Eq. (B.15), the vector 
is written as follows 

^ = e - i 2 rV> 0 . (B.18) 


C. HELLINGER ANGLE AND WYSI: SINGLE QUBIT 
CASE 


(Sf)ji = cos( 2 y)Sji + [1 - cos( 2 y)]( 2 : (0 ) j(Xf)i + sin( 2 y)(A ^) j7 

(B.l 1) 


In this section we provide an explicit calculation of the 
Hellinger angle and Wigner-Yanase skew information for a 
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mixed single qubit state. To achieve these results we will ob¬ 
tain analytical expressions of the inverse matrix, determinant 
and square root for that state. Consider a hermitian contrac¬ 
tive operator 77, i.e., ||77|| < 1, where ||... || defines the opera¬ 
tor (or bound) norm. In this case, the positive operator 7+77 
is invertible and its inverse 

(7 + 77) 1 = 7 - 77 + 77 2 - 77 3 + ... 

= (7 - 77)(7 + 77 2 + 77 4 + ...) 

= (7 - 77)(7 - 77 2 r I (C.19) 

defines a convergent Neumann series [38]. Let us consider 
now a single qubit mixed state p M — (l/2)(7 + r M ■ cr) with 
p e {0, p] and choose FI — r M ■ cr. Here 7 denotes the 2x2 
identity matrix, & - fen, cr?, <+ 3 ) and 7 , is a 3-dimensional 
vector which fulfils Ifyl 2 < 1. By using the vector identity 

. . . 1 , . 

(a ■ cr)(b ■ cr) — (a ■ b)I + i(a x b) ■ <7, is straightforward to 
verify 77 2 = (r M ■ cr) 2 — |iy,| 2 7 and Tr(77) = 0. Given that 

||77|| = III77IH, where |77| = V/77/ = |^|7, in our case follows 
||77|| = ICill|7|| = < 1 and according to Eq. (C.19) we obtain 

P ^ 1 = 2(7 + r M ■ (?)“* = -—7^(7 - y ■ (?) ■ (C.20) 

1 - M 2 

Note that the previous result is singular if the state is a pure 
one. Actually, in this case the inverse operation requires 
another approach known as generalized inverse or Moore- 
Penrose inverse [49]. Returning to the mixed case, it is a sim¬ 
ple task to recognize 1 - l^l 2 as the determinant of the state 
f> t j starting from the identity 

det(7 + 77) = e TllHI+n)] . (C.21) 


According to Eq. (C.20) it can be verified that 

-1 2 _ „ 

(Pn + xl) = (7 + v M • cr) 

2(7 - v u ■ (?) 

" (1 + 2 x)(l - |v^| 2 ) 

_ 2[(1 + 2x)I - y ■ c?] 

(1 + 2x) 2 - |i),[ 2 ’ 


(C.25) 


with )), = (!+ 2x)v /1 , and thus 


PtiiPn + xl) 1 


1 + 2 x - |i >| 2 2x 

(1 + 2x) 2 - + (1 + 2x) 2 - ' ’ 


(C.26) 

Substituting the previous result into Eq. (C.24) and perform¬ 
ing the calculation of both integrals, we finally obtain 


y/pk = ——[c+7 + c-(r„ ■ <?)] 

2V2 


(C.27) 


where 


c* := ^1 + 1^1 ± 7 1 - 1^1 ■ (C.28) 

As pointed out in the main text, the Hellinger angle is deter¬ 
mined by the equation cos[7(po,p^)] = Tr( yfpo yfp^). Start¬ 
ing from the result indicated in Eq. (C.27), we conclude 

cos[X(po.P P )] = ^ [cjc+ + CqC~{ r v ■ r 0 )] . (C.29) 

Analogously, the Wigner-Yanase skew information 
lip^Ky) = —(1 /2)Tr([ -y r p^,K lfi ] 2 ) also depends on the 
square root of the density operator. Considering the local 
observable — m{al + ■ cr) it is possible to prove that 


In fact, since ||77|| < 1 and taking the Taylor series expansion 
ln(l + x) = — 2 *Li(—*)*/& f° r \x\ < 1 , follows 

°° (- 1 )* 

Tr[ln(7 + 77)] = -^--^Tr(77*) 

k=l 


-Z 


" '41" 


k=l 


= ln(l - |?/) , 


(C.22) 


where we used Tr(77 2i+1 ) = 0 and Tr(77 2<: ) = 2|i),| 2i and col¬ 
lected separately even and odd contributions in the infinite 
sum. Therefore, we get 


det(7 + ?„■(?) = !- |r „| 2 . 


(C.23) 


In order to calculate the square root of the density operator 
p^, it is convenient to remember the integral representation 
presented in Eq. (A.2) choosing now A = p M and s = 1/2, i.e., 

yfpii = - f ^pP^(p M + xl)~ l . (C.24) 

n Jo y/x 


and also 


I y/ffy, Ky\ = i —p-(r^ X h v ) ■ & (C.30) 

V2 


[yfp^,^} 2 = -\-{mc ) 2 |r^ x n^| 2 7 , (C.31) 


where we used Eq. (C.27) in order to write 

TxJ 

yJfyK.p = ^=[ac v + c v (f v ■ h tp )]I+ 


m 


2 y/2 


[c^hp + ac^y + ic^iy x h v )] ■ & (C.32) 


and 


y[p; = ^-^[ac + v + c^y • n v )]7+ 


m 

2 y/2 
m 

2 V2 


[c+h 9 + ac^ - iCylfy x h v )] ■ & . (C.33) 


Therefore, the Wigner-Yanase skew-information is given by 
Kp^Ky) = i(wc“) 2 \r v x n ^| 2 . (C.34) 
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Remember that the quantum system evolves under an uni¬ 
tary transformation which does not change the absolute value 
of Bloch sphere radius, i.e., \ff = To = tp. Therefore, 

since c* = cf t and defining i'± - 1 ± J1 - r q, the cosine 
of the Hellinger angle and Wigner-Yanase skew information 
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